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ABSTRACT 

A significant fraction of Lya emitting galaxies (LAEs) at z ^ 5.7 have rest-frame 
equivalent widths (EW) greater than ~ 100A. However only a small fraction of the Lya 
flux produced by a galaxy is transmitted through the IGM, which implies intrinsic Lya 
EWs that are in excess of the maximum allowed for a population-II stellar population 
having a Salpeter mass function. In this paper we study characteristics of the sources 
powering Lya emission in high redshift galaxies. We propose a simple model for Lya 
emitters in which galaxies undergo a burst of very massive star formation that results 
in a large intrinsic EW, followed by a phase of population-II star formation with a 
lower EW. We confront this model with a range of high redshift observations and 
find that the model is able to simultaneously describe the following eight properties 
of the high redshift galaxy population with plausible values for parameters like the 
efficiency and duration of star formation: i-iv) the UV and Lya luminosity functions 
of LAEs at z=5.7 and 6.5, v-vi) the mean and variance of the EW distribution of Lya 
selected galaxies at z=5.7, vii) the EW distribution of i-drop galaxies at z~6, and 
viii) the observed correlation of stellar age with EW. Our modeling suggests that the 
observed anomalously large intrinsic equivalent widths require a burst of very massive 
star formation lasting no more than a few to ten percent of the galaxies star forming 
lifetime. This very massive star formation may indicate the presence of population- 
Ill star formation in a few per cent of i-drop galaxies, and in about half of the Lya 
selected galaxies. 
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1 INTRODUCTION 

Narrow band searches for redshifted Lya lines have 
discovered a large number of Lya emitting galaxies 
with redshifts between z — 4.5 and z = 7.0 (e.g . 
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prominence of the Lya line is quantified by its equivalent 
width (EW), defined as the total flux of the Lya line, F^ ya 
divided by the flux density of the continuum at 1216 A: 
EW= i*Lya//i2i6- Throughout this paper we refer to the 
rest-frame EW of the Lya line (which a factor of (1 + z) 
lower than the EW in the observers frame). 

Approximately 50% of Lya emitters (hereafter LAEs) 
at z = 4.5 and z = 5.7 have lines with EW~ 100 - 500 A 
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J Dawson et al ] |2004l ; lHu et~aU 2004; Shi masaku et al ] |2006h . 
For comparison, theoretical studies conclude that the maxi- 
mum EW which can be produced by a conventional popula- 
tion of stars is 200-300 A. Moreover, this maximum EW can 
only be produced during the first few million years of a star- 
burst, while at l ater times the luminous phase of Lya EW 
gradu ally fades (|Charlot fc Falll 1 19931 ; iMalhotra fc Rhoadsi 
2002). Therefore, observed EWs lie near the upper envelope 
of values allowed by a normal stellar population. 

The quoted value for the upper envelope of EW~ 200 — 
300 A corresponds to the emitted Lya flux. However not 
all Lya photons are transmitted through the IGM, and we 
expect some attenuation. Within the framework of a Cold 
Dark Matter cosmology, gas surrounding galaxies is signifi- 
cantly overdense ; and po ssesses an infall velocity relative to 
the mean IGM (|Barkana|[200i '). As a net result, the IGM 
surrounding high redshift galaxies is significantly opaque to 
Lya photons. Indeed it can be shown that for reasonable 
model assumptions, only ~ 10 — 30% of all Lya phot ons are 
transmitted through the IGM jPiikstra et alj|2007r i. As a 
result, the intrinsic Lya EW emitted by high redshift LAEs 
is systematically larger than observed. Indeed, this observa- 
tion suggests that a significant fraction of LAEs at z ^ 4.5 
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have intrinsic EWs that are much larger than can possibly 
be produced by a conventional population of young stars. 

One possible origin for this large EW population is pro- 
vided by active galactic nuclei (AGN), which can have much 
large r EWs due to their harder spectra (e.g. I Chariot fc Falll 
1993). However, large EW LAEs are not A GN for several 
reaso ns: (1) the Lya lines are too narrow (jDawson et al.l 
[2004|) (2) these objects typically lack high-ionisation state 
UV emission lines, w hich are symptomatic of AGN activity 
ijDawson et ail [2004). and (3) deep X-Ray observations of 
101 Lya emitters by Wang et al. (2004, also see Malhotra et 
al. 2003, Lai et al, 2007) revealed no X-ray emission neither 
from any individual source, nor from their stacked X-Ray 
images. 

Several recent papers have investigated the stel- 
lar content of high-redshift LAEs by comparing stellar 
synthesis models with the observed broad band colors 
i|Finkelstein et alj[2007h . These comparison s are often aided 
by deep IRAC ob servations on Spitzer (|Lai et al.l 120071; 
iPirzkal et a"ill2007h . In this paper we take a different ap- 
proach. Instead of focusing on individual galaxies, our goal 
is to provide a simple model that describes the population 
of Lya emitting galaxies as a whole. This population is de- 
scribed by the rest-frame ultraviolet (UV) and Lya lumi- 
nosity functions (LFs) at z = 5.7 and z = 6.5 , and the 
Lya EW distribution at z = 5.7 (|Shimasaku et all l200fj; 
iKashikawa et al" I l2006h . The sample of high-redshift LAEs 
is becoming large enough that meaningful constraints can 
now be placed on simple models of galaxy formation. 

The outline of this paper is as follows: In § [2] -§ [5] 
we describe our models. In § [S] we discuss our results, 
and compare with results from stellar synthesis models, be- 
fore presenting our conclusions in § [TJ The parameters for 
the background cosmology used throughout this paper are 
Q m = 0.24, = .76, Q b = 0.044, h = 0.73 and cr 8 = 0.74 
jSpergel et alj|2007h . 



2 THE MODEL 

iDiikstra et al] l|2007bl ) found that the observed Lya LFs at 
z = 5.7 and z = 6.5 are well described by a model in which 
the Lya luminosity of a galaxy increases in proportion to 
the mass of its host dark matter Mtot- One can constrain 
quantities related to the star formation efficiency fro m such 
a model (also see lMao et"aill2007l : Tstark et alj|2007f ). 

However, it is also possible to obtain constraints from 
the rest-frame UV-LFs. In contrast to the Lya LF, the UV- 
LF is not affected by attenuation by the IGM, which allows 
for more reliable constraints on quantities related to the star 
formation efficiency. In the first part of this paper (§ [3}§ [p 
we present limited modeling to illustrate parameter depen- 
dences, using the UV-LF to constrain model parameters re- 
lated to star formation efficiency and lifetime. These model 
parameters may then be kept fixed, and the Lya LFs and 
EW distributions used to constrain properties of high red- 
shift LAEs such as their intrinsic Lya EW and the fraction 
of Lya that is transmitted through the IGM. Later, in § O 
we present our most general model, and fit to both the UV 
and Lya LFs, as well as the EW distribution, simultane- 
ously, treating all model parameters as free. 



3 MODELING THE UV AND LYa 
LUMINOSITY FUNCTIONS. 

3.1 Constraints from the UV-LF. 

We begin by presenting a simple mode l for the UV-LF 
|Wvithe fc Loebl 120071 : IStark et ail 12007m . In Figure [TJ we 
show th e rest-frame UV-LFs of LAEs at z = 5.7 and 
z = 6.5 (|Shimasaku et ail 120061 : IKashikawa et abll2006t ). We 
use the following simple prescription to relate the ultravi- 
olet flux density emitted by a galaxy, /1350, to the mass 
of its host dark matter M to t- The total mass of baryons 
within a galaxy is (f2;,/f2 m )Mtot, of which a fraction /» 
is assumed to be converted into stars over a time scale of 
i S ys = eDciimb- Here, edc is the duty cycle and thub(z), the 
Hubble time at redshift z. This prescription yields a star 
formation rate of A/* = f*(yib/Q. m )M /t sys . The star forma- 
tion rate can then be converted into /1350 using t he relation 
/1350 = 7 x 10 27 (M„/[M Q /yr]) erg s _1 Hz" 1 (|Kennicuttl 
1998). The precise relation is uncertain but differs by a fac- 
tor of less than 2 between a normal and a metal-free stel- 
lar population (see e.g. Loeb et al. 2005). Uncertainty in 
this conversion factor does not affect our main conclusions. 
The presence of dust would lower the ratio of /1350 and M„ , 
which could be compe nsated for by increasing /*. However, 
iBouwens et ail (|2006h found that dust in z = 6.0 Lyman 
Break Galaxies (LBGs) attenuates the UV-flux by an aver- 
age factor of only 1.4 (and dust obscuration may be even 
less important in LAEs, see § [6T3J ■ Since this is within the 
uncertainty of the constraint we obtain on /„ , we ignore ex- 
tinction by dust. The number density of LAEs with UV-flux 
densities exceeding /1350 is then given by 

N(> /rase) = enc f°° dM-^-, (1) 
Jm uv ">m 

where Muv is the mass that corresponds to the flux den- 
sity, /1350 (through the relations given above). The func- 
tion dn/dM is the Pre ss-Schechter (1974 ) mass function 
(with the modification of lSheth et aillioOll ) , which gives the 
number density of halos of mass M (in units of comoving 
Mpc~ 3 j3- The free parameters in our model are the duty 
cycle, £dc, of the galaxy, and the fraction of baryons that 
are converted into stars, /„ . We calculated the UV-LF for 
a grid of models in the (edc, /*)-plane, and generated like- 
lihoods C[P] = exp[— 0.5x 2 ], where \ 2 = X^f rdata ( m °deli — 
data;) 2 /of, in which data; and <Ji are the i th UV-LF data 
point and its error, and modeh is the model evaluated at the 
i th luminosity bin. The sum is over iVdata = 8 data points. 
The inset in Figure[TJshows the resulting likelihood contours 



1 Our model effectively states that the star formation rate in a 
galaxy increases linearly with halo mass. This is probably not cor- 
rect. To account for a different mass dependence we could write 
the star formation rate as M* tx M@, where is left as a free pa- 
rameter. However, the range of observed luminosities span only 1 
order of magnitude, and we will show that the choice /3 = 1 pro- 
vides a model that describes the observations well. Furthermore, 
the duty cycle £dc m &y be viewed as the fraction of dark mat- 
ter halos that are currently forming stars. The remaining fraction 
(1 — cue) of halos either have not formed stars yet, or are evolv- 
ing passively. In either case, the contribution of these halos to the 
UV-LF is set to be negligible. 
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Figure 1. Constraints on the star formation efficiency from the 
observed rest-frame UV-luminosity functi ons of LAEs at z = 5.7 
(red squares) and z = 6.5 (blue circles) dShimasaku et al1l2006l ; 
iKashikawa et alj liopr? ). In our best-fit model, a faction /» ~ 
0.06 of all baryons is converted into stars over a time-scale of 
£DC*hub ~ 0-03 Gyr (see text). The inset shows likelihood con- 
tours in the (erjc /* )-plane at 64%, 26% and 10% of the peak 
likelihood. Also shown on the upper horizontal axis is the mass 
corresponding to Mab 1350 m the best-fit model. 



in the (edc, /*)-plane at 64%, 26% and 10% of the peak like- 
lihood. The best fit model has (edc, /*) = (0.03, 0.06) and is 
plotted as the solid line. In the following sections we assume 
this combination of /„ and enc- 



3.2 Constraints from the Lya LF. 

We next model the Lya LF, beginning with the best-fit 
model of the previous section. The number density of LAEs 
at redshi ft z with Lya lumino sities exceeding T a x L a is 
given by (|Diikstra et alj|2007bl ) 



N(> TaX L a ,z) = £DC 



(2) 



where the Lya luminosity and host halo mass, M a are re- 
lated by 



T~ a X L a — £a ' 



i S ys(yr) 



(3) 



In this relation, T a is the IGM transmission multiplied by 
the escape fraction of Lya photons from the galaxy, and 
C a = 2.0 x 10 42 erg s' 1 /(M Q yr" 1 ), is the Lya luminos- 
ity emitted per unit of star formation rate (in Mq yr" 1 ). 
Throughout, C a ,4,2 denotes jC a in units of 10 42 erg s _1 /(Mg 
yr -1 ). We have taken C a ,42 = 2.0, which is appropriate for 
a metallicity of Z = 0.05^0 and a Salpeter IMF (Dijkstra 
et al, 20 07). Note that when comparing to observed l umi- 
nosities |Shimasaku et aLll2006l ; [Kashikawa et al.ll2006l ). we 
have replaced L a with T a x L a . This is because the observed 
luminosities have been derived from the observed fluxes by 
assuming that all Lya emerging from the galaxy was trans- 
mitt ed by the IGM, whe reas there is substantial absorption 
(e.g. lDiikstra et~al ] |2007i ). The product T a x L a may be writ- 
ten as T a x L a = Aixd\(z) S a , where S a is the total Lya flux 
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Figure 2. Joint constraints on the IGM transmission from the 
observed Lya and res t-frame UV LFs of LAEs at z = 5.7 
dShimasaku et al.l2006T) and z = 6.5 [IKashikawa et al.l2006T) . The 
red squares and blue circles represent the data at z = 5.7 and 
z = 6.5. Using the model that best describes the UV-LFs with 
(/star^Dc) = (0.06,0.03) we fit to the observed Lya LF. The 
only free parameter of our model was the fraction of Lya that 
was transmitted through the IGM at z = 5.7, ^,57 (see text). 
Shown in the inset is the likelihood for 7^,57, normalised to a 
peak of unity. The figure shows that in order to simultaneously 
fit the Lya and UV-LFs, only ~ 30% of Lya photons are trans- 
mitted through the IGM. The best fit models are overplotted as 
the solid lines. 



detected on earth and dh(z) is the luminosity distance to 
redshift z. The product T a x L a may therefore be viewed 
as an effective luminosity inf erred at earth. F u rtherm ore, 
the selection criteria u sed by IShimasaku et all l|2006h and 
Kashikaw a et al.l (|2006h limits these surveys to be sensitive 
to LAEs with EW JjlOA. In the reminder of this paper, the 
EW of model LAEs is always larger than this EW m i n , and 
we need not worry about selection effects when comparing 
our model to the data. 

In this section, we set the transmission at z — 6.5 (de- 
noted by 7^,65) to be a factor of ~ 1.2 loweflthan at z — 5.7 
(denoted by 7^.57). Th i s rati o is the median of the range 
found bv lDiikstra et al l (|2007l ). We then calculated the Lya 
LF for a range of T a ^7, and generated likelihoods C[P] = 
exp[— 0.5x 2 ], where \ 2 = J^f^^ (modeL — datai) 2 /cr 2 , for 
each model. Here, data^ and at are the i th data point and its 
error, and model; is the model evaluated at the i th luminos- 
ity bin. The sum is over A^data = 6 points at each redshift. 
In Figure[5]we show the Lya luminosity functions at 2 = 5.7 
and z = 6.5. The red s q uares and blue circl es represent data 
from I Shimasaku et al.1 l|2006l . 2 = 5.7) and IKashikawa et al.l 
(2006, z = 6.5), respectively. 

The likelihood for 7^57 (normalised to a peak of unity) 
is shown in the inset. The best fit model is overplotted 
as the solid lines, for which the value of the transmission 
is 7^,57 = 0.30. The modeling presented in this and the 



2 For our primary results in § \E\ we allow this ratio to be a free 
parameter. The results presented in this section is not sensitive 
to the precise choice of the ratio of IGM transmission at z = 5.7 
and 2 = 6.5. 
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previous section therefore suggests that, in order to simul- 
taneously fit the Lya and UV luminosity functions, only 
~ 30% of the Lya can be transmitted through the IGM. This 
tr ansmission is in good agreement with the results obtained 
bv lDiikstra et afl l|2007h . who modeled the transmission di- 
rectly and found that for reasonable model parameters the 
transmission must lie in the range 0.1 <ST a <J0.3. 

3.3 The Predicted Equivalent Width 

While in agreement with the observed LFs, the model de- 
scribed in § 13.21 does not reproduce the very large ob- 
served equivalent widths. The Lya luminosity can be rewrit- 
ten in terms of the star formation rate (M,) and EW as 
L a = 2.0 x 10 42 erg s _1 M»(M /yr)(EW/16O A). Here we 
have used the relation L a =EWx [u a f(v a )/X a ], where f(v a ) 
is the flux density in erg s _1 Hz -1 at u a , and where we 
denoted the Lya frequency and wavelength by v a and A Q 
respectively. Furthermore, we assumed the spectrum to be 
constant between 1216 A and 1350 A. For L a ^i = 2.0, we 
find a best fit model that predicts LAEs to have an observed 
EW of T a x 160 A ~ 50 A. This value compares unfavor- 
ably with the observed sample that includes EWs exceed- 
ing ~ 100 A in ~ 50% of cases for Lya selected galaxies 
jShimasaku et alj|2006h . 

Thus although our simple model can successfully repro- 
duce the observed Lya and UV luminosity functions, the 
model fails to reproduce the observed large EW LAEs. This 
discrepancy cannot be remedied by changing the intrinsic 
Lya emissivity of a given galaxy, C a '- increasing C a would 
be simply be compensated for by a lower T a and vice versa. 
In the next section we discuss a simple modification of this 
model that aims to alleviate this discrepancy. 



4 THE FLUCTUATING IGM MODEL 

The model described in § [3] assumed that the Lya flux of 
galaxies was subject to uniform attenuation by the IGM. 
In this section we relax this assumption and investigate the 
predicted EWs in a more realistic IGM where transmission 
fluctuates between galaxies. We refer to this model as the 
'fluctuating IGM' model. In this model, a larger transmis- 
sion translates to a larger observed equivalent width. As a 
result, galaxies with large T a are more easily detected, and 
the existence of these galaxies may therefore affect the ob- 
served EW-distribution for Lya selected galaxies, even in 
cases where they comprise only a small fraction of the in- 
trinsic population. In this section we investigate whether 
this bias could explain the anomalously large observed EW. 

We assume a log-normal distribution for 7^,57, 

P{u)du = x exp ( ~ (M ~ <M>)2 W (4) 

where (it) = log(7^,57) is the log (base 10) of the mean 
transmission and a u is the standard deviation in log-space. 
Throughout this section we drop the subscript '57'. Eq ((4]) 
may be rewritten in the form 



/(> u) = 




1 — t 




Figure 3. Same as Figure 2. However, instead of assuming a 
single value of IGM transmission 7^,57, we assumed a log-normal 
distribution of IGM transmission with a mean of log 7^57 and 
standard deviation of (in the log) a u (see text). This reflects the 
possibility that the IGM transmission fluctuates between galaxies. 
The inset shows likelihood contours for (log T a .57, c u ). Increasing 
a u flattens the luminosity function (and moves it upward), which 
is illustrated by the model LFs shown as solid lines, for which 
we used (7^57, cr u ) = (0.27,0.2) (shown as the thick black dot 
in the inset). The best-fit model to the data has cr u ~ (which 
corresponds to the model shown in Figure 2). 

which gives the fraction of LAEs with log T a > u. The num- 
ber density of LAEs is then given by 

N{>T a xL a ) =e DC I" ' dM-^ 7 f{>u(T a xL a ,M)) (6) 
Jo dM 

where 

u(T a xL a ,M) = log (^y^)- (7) 

Eq ((6} differs from Eq ((2)1 in two ways: (1) there is no lower 
integration limit, and (2) there is an additional term /(> 
u). These two differences reflect the facts that all masses 
contribute to the number density of LAEs brighter than T a x 
L a , and that lower mass systems require larger transmissions 
(Eq[7J) which are less common (Eq[5|. In the limit a u — > 0, 
the function /(> u) 'jumps' from to 1 at M m i n (Eq[3]), 
which corresponds to the original Eq ([2]). 

In this formalism we may also write the number density 
of LAEs with transmission in the range u±du/2 = log T a ± 
f " , which is given by 

N{u)du = P{u)du f dM -JU- ( 8 ) 

Here M m in(w) is the minimum mass of galaxies that can 
be detected with a transmission in the range u ± du/2 (for 
Aftot < Af m in the total flux falls below the detection thresh- 
old). The number density of LAEs with transmission in the 
range u ± du/2 may be used to find the number density of 
LAEs with equivalent widths in the range log EW± dlog 2 EW 
via the relation EW= 1607^ A (for the choice £^,42 = 2.0, 
see § 13 - 2 p . Eq © shows that the observed equivalent width 
distribution takes the shape of the original transmission dis- 
tribution, modulo a boost which increases towards larger 
EWs. 

As in § 13.21 we assume the best-fit model parameters 
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for £dc and /* derived from the UV-LFs determined in 
§ 13.11 We calculate model Lya LFs on a grid of models in 
the (a u , (7^})-plane, and generate likelihoods following the 
procedure outlined in § 13.21 The results of this calculation 
are shown in Figure [3] where we plot the Lya LFs together 
with likelihood contours in the (o~ u , (7^})-plane (inset). The 
best-fit models favor no scatter in T a (a u ~ 0). The reason 
for this is that for any given model, a scatter in T a serves 
to flatten the model LF. However, the observed Lya LF is 
quite steep, and as a result the data prefers a model with no 
scatter. Furthermore, a scatter in T a results in a model LF 
that lies above the original (constant transmission) LF at 
all T a x L a . This also explains the shape of the contours in 
the (<T U , (7^))-plane; increasing a u must be compensated for 
by lowering {%,)). To illustrate the impact of a fluctuating 
T a , the model LFs shown in Figure [3] are not those of the 
best-fit model, but of a model with {(%,), <r u ) = (0.27,0.2). 

Figure [3] shows the observed EW-distribution (solid 
line) at z = 5.7 associated with the Lya LFs shown in Fig- 
ure [3] This model may be compared to the observed distri- 
bution (shown as the histogram, data from Shimasaku et al, 
2006). The upper horizontal axis shows the transmission cor- 
responding to each equivalent width. The dotted line shows 
the distribution of transmission (given by Eq[3J|. Note that 
the range of transmissions shown in Figure [4] extends to 
T a > 1, which of course is not physical. The model EW- 
distribution has a tail towards larger EWs. However, the 
model distribution peaks at EW~ 50 A. As was seen in 
the constant transmission model, this is clearly inconsistent 
with the observa t ions, which favor values of EW~ 100 A. 
Shimasaku et al. I (|2006l ) also show a probability distribution 
of EW, which is probably closer to the actual distribution. 
Although this distribution does peak closer to EW~ 50 A, 
it is still significantly broader than that of the model (see 
§[6]for more arguments against the fluctuating IGM model). 
We point out that the model distribution shown in Figure 2] 
is independent of the assumed value of C a ,42- 



5 LYa EMITTERS POWERED BY VERY 
MASSIVE STARS. 

In § [3] we demonstrated that a simple model where L a and 
Luv were linearly related to halo mass can reproduce the 
UV and Lya LFs, but not the observed EW distribution. 
In § [4] we showed that this situation is not remedied by a 
variable IGM transmission, and that favored models have a 
constant transmission. In this section we discuss an alternate 
model, which leads to consistency with both the observed 
Lya LFs, UV-LFs, and the EW-distribution. In this model, 
galaxies are assumed to have a bright Lya phase (hereafter 
the 'population Ill'-phase) which lasts a fraction fm of the 
galaxies' life-time. After this the galaxy's Lya luminosity 
drops to the 'normal' value for population II star formation. 

This model may be viewed as an extension of th e idea 
originally described by M alhotra fc Rhoadsl ((20 02). that 
large EW LAEs are young galaxies in the early stages of 
their lives. In this picture, the sudden drop in Lya luminos- 
ity could represent i) a sudden drop in the ionising lumi- 
nosity when the first O-stars died, or ii) an enhanced dust- 
opacity after enrichment by the first type-II supernovae. Al- 
ternatively, our parametrisation could represent a scenario 
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Figure 4. Comparison of the observed equivalent width distribu- 
tion (EW, histogram), with the model prediction for a model in 
which we assumed a log-normal distribution of IGM transmission 
with (7^,57, o- u ) = (0.27, 0.2) (see Fig 3). The EW is related to T a 
via EW= 1607^. A. The dotted line shows the fraction, /(> T a ) 
(shown on the right vertical axis), of galaxies with a transmission 
greater than T a (EqE). Galaxies with large T a are more easily 
detected, hence the large T a (EW) end is boosted considerably, 
resulting in closer agreement (but not close enough) to the data. 

in which the population III phase ended after the first pop- 
ulation III stars enriched the surrounding interstellar gas 
from which subsequent generations of stars formed. Hence, 
we refer to this model as the 'population III' model. We will 
show that to be consistent with the large values of the ob- 
served EW, a very massive population of stars is required 
during the early stages of star formation. 

To minimise the number of free parameters we modeled 
the time dependence of the Lya EW as a step-function. The 
number density of LAEs is then given by 

f°° dn 

N{> T a x L a ,z) = fm x e DC / dM ^jj( z ) + ( 9 ) 

J M a IU dm 

(l-/„:)xe DC £^M-^( Z ). 

Here, M a ,n is the mass related to T a x L a through T a xL a = 
C a xT a x M* , while M a ,m is the population III mass, which 
is calculated with C a replaced by L a = (EWm/160 A) x 2 x 
10 42 erg s _1 . 

Whereas in previous sections we chose fiducial or best 
fit parameters for illustration, for the model described in this 
section we take the most general approach. We fit the model 
simultaneously to the UV-LF and Lya LFs, as well as to 
the observed EW-distribution of Lya selected galaxies. This 
model predicts two observed equivalent widths (T a x EWni 
and T^xEWn) in various abundances. The associated mean 
and variance from the model are compared to the observed 
EW-distribution, which has a mean of (EW) = 120 ± 25 A, 
and a standard deviation of ctew = 50 ± 10 A. 

Our model has 6 free parameters 
(eDC,/»,7" Q ,57,T Q ,65,/ni,EWni). We produce likelihoods 
for each parameter by marginalising over the others in 
this space. The lower set of panels in Figure [5] show 
likelihood contours for our model parameters at 64%, 
26% and 10% of the peak likelihood. The best-fit models 
have EWm ~ 600 - 800 A and f m = 0.04 - 0.1 which 
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Figure 5. Marginalised constraints on the 6 population-Ill model parameters, (cdc, /* ,T a ,57, T a ,65, fill, EWm), are shown in the lower 
three panels. In the best-fit population-Ill model, each galaxy goes through a luminous Lya phase during which the equivalent width 
is EWin = 600 — 800 A for a fraction fm = 0.04 — 0.1 of the galaxies life time. Although the bright phase only lasts a few to ten per 
cent of their life-time, galaxies in the bright phase are more easily detectable, and the number of galaxies detected in Lya surveys in the 
bright phase is equal to that detected in the faint phase. This is also demonstrated by the thick solid line (with label '1.0') in the lower 
left panel, which shows the combination of fm and EWni that produces equal numbers of galaxies in the population III and II phase 
(the dashed lines are defined similarly, also see Fig [6jl . The best-fit intrinsic equivalent widths are in excess of the maximum allowed for 
a population-II stellar population having a Salpeter mass function. Therefore, this model requires a burst of very massive star formation 
lasting no more than a few to ten percent of the galaxies star forming lifetime, and may indicate the presence of population-Ill star 
formation in a large number of high-redshift LAEs. 



corresponds to a physical timescale for the population-Ill 
phase of fm x e D c X ihub ~ 4 — 50 Myr (for 0.1 <CeDC <^0.5). 
The model Lya luminosity functions at z=5.7 and 
z=6.5 described by (edc, /*, T a ,57, T a ,a5, fm, EWm) = 
(0.2,0.14,0.22,0.19,0.08,650 A) are shown as solid lines 
and provide good fits to the data. The model produces 
two observed EWs, namely Z^xEWn = 35 A and T a x 
EWm = 143 A . It is worth emphasising that the emitted 
EW of the bright phase depends on the choice C a ,42 via 
EWm = 650(£ a ,42/2.0) A. Hence, a lower/higher value 
of C a ,42 would decrease/increase the intrinsic brightness 
of the 'population IIP phase. Note that €. a ,42 = 1.0 when 
LAEs formed out of gas of solar metallicity, which is 
unreasonable given the universe was only ~ 1 Gyr old at 
2 = 6. Furthermore, £ a ,42 = 1 would have yielded a best-fit 
7^,57 = 0.5, which i s well outside the range calculated by 
Dijks tra et al.l (|2007l ). We there conclude that C a ,42 is m 
excess of unity. 

In performing fits we have fixed the value of EWn to 
correspond to a standard stellar population, and then ex- 
plored the possibility that there might be a second phase of 
SF producing a larger EW. Our modeling finds strong sta- 
tistical evidence for this early phase and rules out the null- 
hypothesis that properties can be described by population-II 



stars alone at high confidence (grey region in the lower left 
panel inset of Fig [5}. Despite the fact that the best fit model 
has a bright phase which lasts only a few per cent of the total 
star formation lifetime, the two populations of LAEs are sim- 
ilarly abundant in model realisations of the observed sample 
in Lya selected galaxies (see § 15.11 for a more detailed com- 
parison to the observed EW distributions). This is shown 
in the lower left panel in Figure [5] in which the solid line 
(with label '1.0') shows the combination of fm and EWm 
for which the observed number of galaxies in the popula- 
tion III phase (Mii) equals that in the population II phase 
(Nix). The dashed lines show the cases Nui/Nu = 0.3 and 
Niii/Nu = 3.0. The duration of the bright Lya phase meets 
theoretical expectations for a burst of star formation, while 
the large EW requires a very massive stellar population (e.g 
ISchaererl 2003 ; iTumlinson et al.ll2003l ) . In summary, in order 
to reproduce both the UV and Lya LFs, and the observed 
population of large EW galaxies, we require a burst of very 
massive star formation lasting <T0 per cent of the galaxies 
lifetime. 
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5.1 EW Distribution of UV and Lya Selected 
Galaxies 

IStanwav et all ()2007l ) show that 11 out of 14 LAE candidates 
among i-drop galaxies in the Hubble Ultra Deep Field have 
EW< 100 A. If galaxies are included for which only upper or 
lower limits on the EW are available, then this fraction be- 
comes 21 out of 26. Thus the distribution of EWs for i-drop 
selected gal axies differs strongly fro m the EW-distribution 
observed bv lShimasaku et alj 1 20061 ). We next describe why 
this strong dependence of the observed Lya EW distribu- 
tion on the precise galaxy selection criteria arises naturally 
in our population III model. 

We first assume that the Lya selected and UV-selected 
galaxies were drawn from the same population (this assump- 
tion is discussed further in § I6.3[) . In our model a galaxy that 
is selected based on it's rest-frame UV-continuum emission 
has a probability fm of being observed in the Lya bright 
phase, while the probability of finding a galaxies in the 
Lya faint phase is 1 — fm. In § [5] we found fm ~ 0.1, 
hence an i-drop galaxy is ~ 10 times more likely to have 
a low than a high observed EW. If we denote the number 
of galaxies with EW> 100 A by Nm, and the number of 
galaxies with EW< 100 A by Nit, then the model predicts 
Nm/Nu =/iii/(1 — /in) ~ 0.1, while the observed fraction 
including the galaxies for which the EW is known as upper 
or lower limit is Nm/Nu = 0.19 ±0.05. Therefore the quali- 
tative difference in observed Lya EW distribution among i- 
drop galaxies in the HUDF and among Lya selected galaxies 
follows naturally from our two-phase star formation model. 
Note that our model predicts population III star formation 
to be observed in /m/(l - fm) ~ 10% of the z = 6.0 LBG 
population. 

The dependence of the observed EW distribution on the 
selection criteria used to construct the sample of galaxies is 
illustrated in Figure [6] To construct this figure, we have 
taken the best-fit population III model of § [5] For the pur- 
pose of presentation, we let the IGM fluctuate according to 
the prescription of §[4] with a u — 0.1, so that the model pre- 
dicts a finite range of EWs in each phase. The left and right 
panels show the predicted EW distribution for Lya selected 
(left panel) and UV-selected (right panel) galaxies as the 
solid lines, respectively. For a UV-selected galaxy the prob- 
ability of being in the bright phase and having an observed 
EW in the range EWm x (T a ±dT a /2) is fmP(T a )dT a . Here 
P(T a )dT a is the probability that the IGM transmission is in 
the range T a ± dT a /2, which is derived from Eq[3] The units 
on the vertical axis are arbitrary, and chosen to illustrate 
the different predicted and observed Lya EW distribtions 
for the two samples at large EWs. The observed distribu- 
tions for Lya a nd UV selected galaxies, show n as histograms, 
are ta ken from Shimasa ku et al.l (|2006l ) and IStanwav et al.l 
l|2007l ). respectively. Figure [5] clearly shows that both the 
predicted and observed Lya selected samples contain signif- 
icantly more large EW LAEs than the UV-selected sample. 
Our model naturally explains the qualitative shape of these 
distributions and their differences. 

Before proceeding we mention a caveat to the distribu- 
tions shown in Figure [6] In our model all galaxies have an 
EW of T a xEWn ~ 25—35 A during th e population II phase, 
while in contrast IStanwav et alj (|2007h do not detect 10 out 
of 26 LBGs, which implies that ~ 40% of LBGs have an 




T a x EW T a x EW 



Figure 6. Comparison of the predicted EW distribution for UV 
and Lya selected galaxies. The best-fit population-Ill model (see 
§ [5j was used. In order to get a finite range of observed EWs (in- 
stead of only two values at T a xEWn and T a xEWni, we assumed 
the IGM transmission to fluctuate. The units on the vertical axis 
are arbitrary. The figure shows that in our population III model, 
the Lya selected sample contains a larger relative fraction of large 
EW LAEs than the UV-selected (i-drop) sample, which is quali- 
tatively in good agreement with the observations (shown by the 
histograms) . 

EW <J6 A. Thus there is a descrepancy between our model 
and the observations with respect to the value of EW in the 
population II phase. The resolution of this discrepancy lies 
in the fact that the very low EW emitters are drawn from 
the UV (i-drop) sample and not the Lya selected sample 
our model was set up to describe. This issue is discussed in 
more detail in § 16.31 

An E W distribution of dropout sources was also pre- 
sented by iDow-Hygelund et al.l (120071). These a uthors per- 
formed an analysis similar to IStanwav et al.l (|2007T) and 
found 1 LAE with EW= 150 A among 22 candidate 
z=6.0 LBGs. When interpreted in reference to our model, 
this translates to Nm/Nu ~ 5%, which is consistent 
with the model predictions. Therefore, when interpreted 
in light of a two-phase star formation history and differ- 
ent selection methods, the EW distribution observed by 
IDow-H ygelund ct alj (120071 ) is consistent with that found 
bv lShimasaku et all (|2006l ). 

If population III star formation does provide the ex- 
planation for the very large EW Lya emitters, then we 
would expect the large EW emitters to become less common 
with time as the mean metallicity of the Universe increased. 
To test this idea, we can compare the EW- distribution at 
z = 5. 7 with the results at lower redshift from lShaplev et ail 
l|2003l ) who found that ^0.5% of z = 3 LBGs have Lya 
EW> 150 A , and that <;2% of z = 3 LBGs to have Lya 
EW> 100 A. IDow-Hygelund et alj (|2007n argue that the 
fraction of large EW Ly a lines at z = 6 is co nsistent with 
that observed at z = 3 (|Shaplev et alj [20031 ). However, if 
the EW distribution did not evolve with redshift, then the 
probability that a sample of 22 LBGs will contain at least 
1 LAE with EW ^150 A is <10%. Thus the hypothesis 
that the observed EW distribution remains constant is ruled 
out at t he ~ 90% leve l . On the other hand, in a similar- 
analysis IStanwav et al. I d2007h found 5 out of 26 LBGs to 
have an EW J>100 A. If the EW distribution did not evolve 
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with redshift, then the probability of finding 5 EW J>100 
A in t his sample is only ~ 10~ 4 . Furthermore. iNagao et all 
(j2007l l recently found at least 5 LAEs with EW> 100 A at 
6.0 <.z <:6.5, and conclude that 8% of i'-drop galaxies in 
the Subaru Deep Field have EW> 100 A, which is signifi- 
cantly larger than the fraction of large EW LBGs at z = 3. 
Therefore, the observed EW distribution of LBGs at z — 6 is 
skewed more toward large EWs than at z — 3. The strength 
of this result is increased by the fact that the IGM is more 
opaque to Lya photons at z = 6 than at z = 3. Thus we con- 
clude that the intrinsic EW distribution must have evolved 
with redshift. 



6 DISCUSSION 

6.1 Comparison with Population Synthesis 
Models 

Population synthesis models have suggested that the broad 
band colors of obser ved LAEs are best described with young 
stellar populations dGawiser et af 120061; Pirzkal et al. 2007: 



IFinkelstein et all l2007f ). lLai et all l|2007r ) found thestellar 
populations in three LAEs to be 5 — 100 Myr old, and possi- 
bly as old as 700 Myr (where the precise age upper limit de- 
pends on the assumed star formation hist o ry of t he galaxies). 
However as was arg ued bv lPirzkal et~aH (| 20071 ). since these 
galaxies were selected based on their detection in IRAC, 
a selection bias towards o lder stellar pop ulations may ex- 
ist (al so see lLai et al.ll2007h . Furthermore. IFinkelstein et al.l 
l|2007h found that LAEs with EW> 110 A have ages <;4 
Myr, while LAEs with EW< 40 A have ages between 20- 
400 Myr. This latter result in particular agrees well with 
our population III model. On the other hand, in a fluctu- 
ating IGM model for example, the EW of LAEs should be 
uncorrelated with age. 

In models presented in this paper, on average /* ~ 0.15 
of all baryons are converted into stars within halos of mass 
Mtot ~ 10 10 - lO 11 ^^, yielding stellar masses in the range 
M» = 10 8 -10 9 M Q (|Diikstra et al.ll2007bl ; lst"ark et alj|2007h . 
This compares unfavorably with the typical stellar masses 
found observationally in LAEs which can be as low as 
M* = 10 6 - 1O 7 M (IFinkelstein et all 120071 ; iPirzkal et ail 
120071 ). However, the lowest stellar masses are found (natu- 
rally) for the younger galaxies. Indeed, the LAEs with the 
oldest stellar populations can have stellar masses as large 
a.s 1O 1O M . Thus, we do not find the derived stellar masses 
in LAEs to be at odds with the results of this paper. If 
significant very massive (or population III) star formation 
indeed occurred in high redshift LAEs, then one may ex- 
pect these stars to reveal thems elves in un usual broad-band 
colors (e.g. IStanwav et alj|2005h . However. iTumlinson et all 
l|2003l ) have shown that the most distinctive feature in the 
spectrum of population H I stars is the number of H and He 
ionising photons (also see iBromm et af]|200ll ). Since these 
are (mostly) absorbed in the IGM, the broad band spec- 
trum of population III stars is in practic e difficult to dis- 
tingu ish from a normal stellar population l|Tumlinson et all 
I2003T I. especi ally when nebular conti nuum emission is taken 
into account l|Schaerer fc Polio 2005], see their Fig 1). Hence, 
population III stars would not necessarily be accompanied 
by unusually blue broad band colors. 



6.2 Alternative Explanations for Large EW LAEs 

We have shown that a simple model in which high-redshift 
galaxies go through a population-Ill phase lasting <;15 
Myr can simultaneous ly explain the observe d Lya LFs at 
z = 5.7 and z = 6.5 CKas hikawa et af]|2006l ). and the ob- 
served EW-distribution of Lya selected galaxies at z — 5.7 
(Shimasaku et al. 200(|). In addition, this model predicts the 
much lower EWs found in the population of UV selected 
galaxies (Stanway et al, 2007, see § 15. ip . Moreover the con- 
straints on the population-Ill model parameters such as the 
duration and the equivalent width of the bright phase are 
physically plausible, and consistent with existing population 
synthesis work (see § 16.11) . 

Are there other interpretations of the large observed 
EWs? One possibility was discussed in § [4] where we showed 
that the simple model in which the IGM transmission fluc- 
tuates between galaxies reproduces the LFs, but fails to 
simultaneously reproduce the Lya LFs and the observed 
EW-distribution. I n addition, this mode l fails to reproduce 
other observations. iDiikstra et al.l (|2007h calculated the im- 
pact of the high-redshift reionised IGM on Lya emission 
lines and found the range of plausible transmissions to lie 
in the range 0.1 < T a < 0.3. This work showed that it is 
possible to boost the transmission to (much) larger values 
but not without increasing the observed width of the Lya 
line. Absorption in the IGM typically erases all flux blue- 
ward of the Lya resonance, and when infall is accounted 
for, part of the Lya redward of the Lya resonance as well. 
This implies that Lya lines that are affected by absorption 
in the IGM are systematically narrower than they would 
have been if no absorption in the IGM had taken place. It 
follows that in the fluctuating IGM model, Lya EW should 
be strongly correlated with the observed Lya line width (or 
FHWM). This correlation is not observed. In fact, observa- 
tions suggest that an anti-correlatio n exists between EW 
and FWHM l|Shimasaku et all l200rj ; iTapken et aTll2007l ). 
This anti-correlation provides strong evidence against the 
anomalously large EWs being produced by a fluctuating 
IGM transmission. 

A second possibility is the presence of galaxies with 
strong superwinds. The models of IDiikstra et al.l l|2007h did 
not study the impact of superwinds on the Lya line pro- 
file. The presence of superwinds can cause the Lya line to 
emerge with a systematic redshift relative to the Lya reso- 
nance through back scattering of Lya photons off the far side 
of the shell that surrounds the galaxy (|Ahn et al.1l2003l ; lAhnl 
12004 lHansen fc 6hll2006l ; IVerhamme et alj|2006l ) . However, 
superwinds tend not only to redshift the Lya line, they also 
make the Lya line appear broader than when this scatter- 
ing does not occur. As in the case of the fluctuating IGM 
model, this results in a predicted correlation between EW 
and FWHM, which is not observed. Furthermore in wind- 
models, the overall redshift of the Lya line, and thus T a , 
increases with wind velocity, v w . This predicts that EW in- 
creases with wind velocity. H owever, observations of z — 3 
LBG s bvlShaplev et al.ld2003j ) show that EW correlates with 
v^ 1 l|Ferrara fc Ricotta 120061 ) . We therefore conclude that 
the large EW in LAEs cannot be produced by superwind 
galaxies. 

A third possibility might be that within the Lya emit- 
ting galaxy, cold, dusty clouds lie embedded in a hot inter- 
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Figure 7. Comparison of the best-fit population III model 
(shown in F i gure [5} with the UV-LF constructed by 
iBouwens et al.l hOOdt) using ~ 300 LBGs in the Hubble Deep 
Fields. This good agreement is likely a coincidence since in our 
model all LBGs have an observed EW J>30 A, while the obser- 
vations show that ~ 40% of all LBGs are not detected in Lya 
(EW ^6 A). This overlap could possibly be (partly) due to a 
lower dust content of LAEs relative to their Lya quiet counter- 
parts (see text). 

cloud medium of negligible Lya opacity. Under such condi- 
tions, the continuum photons can suffer more attenuation 
than Lya photons which bounce from cloud to cloud and 
mainly propagate through the hot, transpa rent inter-cloud 
medium (|Neufeldlll99ll ; Irlansen fc OhlbOOrj ). This attenua- 
tion of continuum leads to a large EW. We point out that in 
this scenario, large EW LAEs are not intrinsically brighter 
in Lya. At fixed Lya flux, one is therefore equally likely 
to detect a low EW LAE. In other words, to produce the 
observed EW distribution one requires preferential destruc- 
tion of continuum flux by dust in ~ 50% of the galaxies. 
Currently there is no evidence that this mechanism is at 
work even in one galaxy. Furthermore, the rest-frame UV 
colors of galaxies in the Hubble Ultra Deep Field imply that 
dust in high-redshift galax ies suppresses the con tinuum flux 
by only a factor of ~ 1.4 jBouwens et al.l [20061 ). The maxi- 
mum boost of the EW in a multi-phase ISM is therefore 1.4, 
which is not nearly enough to produce intrinsic equivalent 
widths of EW~ 600 - 800 A. In summary, the only model 
able to simultaneously explain all observations calls for a 
short burst of very massive star formation. 

6.3 Comparison with the LBG Population 

In § 15.11 we have shown that the observed EW distributions 
of Lya selected and i-drop galaxies and their differences can 
be reproduced qualitatively with our population III model. 
However, in our model all high-redshift galaxies have an ob- 
served EW of at least T a x EWn ~ 30 A, whereas many 
i-drop galaxies are not detected in Lya. iKashikawa et all 
(2006) show that the UV-LFs of LAEs a t z = 6.5 and 
z = 5 .7 overlap with that constructed by IBouwens et alj 
(2006) from a sample of ~ 300 z = 6 LBGs discovered in 
the Hubble Deep Fields. Naively, this overlap implies that 



LBGs and LAEs are the same population and therefore that 
all LBGs should be detected by Lya surveys. Since Lya sur- 
veys only detect galaxies with EW ^>20 A, this suggests all 
LBGs should have a Lya EW ;>20 A contrary to observa- 
tion. To illustrate this point further, we have taken the best- 
fit population III model shown in Figure[S]and compared the 
model predicti o ns fo r the rest-frame UV-LF with that of 
IBouwens et alj l|2006l ) in Figure [7] Clearly, o ur best-fit pop- 
ulatio n III model fits the data well. However. IStanwav et al.l 
|2007t ) found ~ 40% of i-drop galaxies in the HUDF to have 
an observed EW <;6 A, and a similar result was presented 
bv lDow-Hygelundet alj i|2007l ). 

Two effects may help reconcile these two apparently 
conflic ting sets of observations: (i) iDow-Hygelund et al.l 
found Lya emitting LBGs to be systematically 
smaller. That is, for a fixed angular size, the zsso-band flux 
of LAEs is systematically higher with Azgso ~ — 1- If we 
assume that the angular scale of a galaxy is determined by 
the mass of its host halo, then this implies that for a fixed 
mass the zsso-band flux of LAEs is systematically higher, 
and (ii) only a fraction of LBGS are LAEs. The drop-out 
technique used to select high redshift galaxies is known to 
introduce a bias against strong LAEs, as a strong Lya line 
can affect the broad band colors of high-redshift galaxies. 
This may cause ~ 10 — 46% of large EW LAEs to be m issed 
using the i-drop technique (|Dow-Hygelund et alj|2007r ). 

If only a fraction f a of all LBGs are detected in Lya, 
then effect (i) would explain why the UV-LF of LAEs lies 
less than a factor of l/f a below the observed UV-LF of 
the general population of LBGs. This is because a more 
abundant lower mass halo is required to produce the same 
UV-flux in LAEs, which would shift the LF upwards. In 
addition, effect (ii) may reduce this difference even further. 
It follows that these two effects combined may cause the LFs 
to overlap. Thus the overlap of the UV and Lya selected UV- 
LFs appears to be a coincidence, and not evidence of their 
being the same population of galaxies. This implies that our 
model is valid for Lya selected galaxies, but not the high- 
redshift population as a whole and explains the lack of very 
low EWs in UV selected samples discussed in § 15.11 

The reason why LAEs may be brighter in the UV for 
a fixed halo mass is un c lear. I t is possibly related to dust 
content. IBouwens et alj |2006l) found the average amount 
of UV exctinction to be 0.4 mag in the total sample of 
z — 6 LBGs. This value is close to the average excess zsso- 
band flux detected from L AEs for a given angular scale 
l|Dow-Hvgelund et ail 120071 ). If LAEs contain less (or no) 
dust, then this would explain why they are brighter in the 
UV and thus why they appear more compact. The possi- 
bility that 'Lya quiet' LBG contain more dust than their 
Lya emitting counterparts is not very surprising, as a low 
dust abundance has the potential to eliminate the Lya line. 
Thus LAEs could be high redshift galaxies with a lower dust 
con tent. 

IShimasaku et all |2006l) and lAndo et all (|2006l ) found 
that luminous LBGs, Mw <; — 21.0, typically do not con- 
tain large EW Lya emission lines. This deficiency of large 
EW LAEs among UV-bright sources is not expected in 
our model, and may reflect that UV-bright sources are 
more massive, mature, galaxies that cannot go through a 
population-Ill phase anymore. It should be pointed out 
though that the absence of large EW LAEs among galax- 
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ies wi th Muv <, — 21.0 in the survey of IShimasaku et al.l 
(2006) is consistent with our model: The observed num- 
ber density of sources with Muv <J — 21.0 is ~ 5 x 10 -5 
cMpc" 3 (see Fig[T}. In our best-fit pop-Ill model (shown in 
Fig©? a fraction /m ~ 0.08 of these galaxies would be in the 
bright phase. This translates to a number density of large 
EW LAEs of ~ 4 x 10 -6 cMpc" 3 . Given the survey volume 
of ~ 2 x 10 5 cMpc 3 , the expected number of large EW LAEs 
with Muv <; — 21.0 is ~ 0.8, and the absence of large EW 
LAE among UV bright sources is thus not surprising. 



6.4 Clustering Properties of the LAEs 

In our model large EW LAEs are less massive by a factor 
of EWiii/EWn ~ 4 at fixed Lya luminosity. Since clus- 
tering of dark matter halos increases with mass, it follows 
that our model predicts large EW LAEs to be clustered less 
than their low EW counterparts (at a fixed Lya luminos- 
ity) . The clustering of LAEs is typically quantified by their 
angular correlation function (ACF), w(8), which gives the 
excess (over random) probability of finding a pair of LAEs 
separated by an angle 8 on the sky. The ACF depends on 
the square of the bias parameter (w(8) oc b 2 (m)), which for 
galaxies in the population II phase is ~ 1.24 — 1.4 times 
larger than for galaxies in the population III phase, for the 
mass range of interest. This implies that the clustering of low 
EW LAEs at fixed Lya luminosity is enhanced by a factor 
of ~ 1.5 — 2.0. Exist i ng det ermi nations of the ACF of L AEs 
bv lShimasaku et all (|2006h and iKashikawa et all (|2006h are 
still too uncertain to test this prediction. 



7 CONCLUSIONS 

Observations of high redshift Lya emitting galaxies (LAEs) 
have shown the typical equivalent width (EW) of the Lya 
line to increase dramatically with redshift, with a signifi- 
cant fraction of the galaxies lying a t z ^ 5.7 hav i ng an 
EW ^100 A. Recent calculations by iDiikstra et al.l (|2007l ) 
show that the IGM at z ^ 4.5 transmits only 10 — 30% 
of the Lya photons emitted by galaxies. In this paper we 
have investigated the transmission using a model that re- 
produces the observed Lya and UV LFs. This model re- 
sults in an empirically determined transmission of T a ~ 
0.30(£ Q ,42/2.0) _1 , where C a ,^2 denotes the Lya luminos- 
ity per unit star formation rate (in Mq yr _1 ) C a in units 
of 10 42 erg s _1 (§[3]). This value is in good agreement with 
earlier theoretical results. 

If only ~ 30% of all Lya that was emitted by high 
redshift galaxies reaches the observer, then this impies that 
the intrinsic EWs are systematically (much) larger than ob- 
served in many cases. To investigate the origin of these very 
high EWs, we have developed semi-analytic models for the 
Lya and UV luminosity functions and the distribution of 
equivalent widths. In this model Lya emitters undergo a 
burst of very massive star formation that results in a large 
intrinsic EW, followed by a phase of population-II star for- 
mation that produces a lower EW(f|. This model is referred 



3 Technically, the model discussed in §[5]only specifies that galax- 
ies go through a 'population-Ill phase for a fraction /ni ~ 0.1 of 



to as the 'population III model' and is an extension of the 
idea originally described by lMalhotra fc Rhoadsi (|2002l ) , who 
proposed large EW Lya emitters to be young galaxies. 

The population III model in which the Lya equivalent 
width is EWm ~ 650(£ Q ,42/2.0) A for ^50 Myr, is able 
to simultaneously describe the following eight properties of 
the high redshift galaxy population: i-iv) the UV and Lya 
luminosity functions of LAEs at z=5.7 and 6.5, v-vi) the 
mean and variance of the EW distribution of Lya selected 
galaxies at z=5.7, vii) the EW distribution of UV- selected 
galaxies at z~6 (§EJ|, an d viii) the observed correlation of 
stellar age and mass with EW (§ 16. lj) . Our modeling sug- 
gests that the anomalously large intrinsic equivalent widths 
observed in about half of the high redshift Lya emitters 
require a burst of very massive star formation lasting no 
more than a few to ten percent of the galaxies star forming 
lifetime. This very massive star formation may indicate the 
presence of population-Ill star formation in a large number 
of high-redshift LAEs. The model parameters for the best-fit 
model are physically plausible where not previously known 
(e.g. those related to the efficiency and duration of star for- 
mation), and agree with estimates where those have been 
calculated directly (e.g the IGM transmission, EWm, and 
/m). 

In addition, we argued that the observed overlap of the 
UV-LFs of LAEs with that of z~ 6 LBGs appears to be at 
odds the observed Lya detection rate in high-redshift LBGs, 
suggesting that LAEs and LBGs are not the same popula- 
tion. A lower dust content of LAEs relative to their 'Lya 
quiet' counterparts would partly remedy this discrepancy, 
and could also explain why LAEs appear to be typically 
more compact (§ 16. 3p . 

Semi-analytic modeling of the coupled reionisation and 
star formation histories of the universe suggests that popu- 
lation III star formation could stil l occur after the bulk of 
reionisation had been completed (IScannapieco et al. I l2003l : 
ISchneider etaH 120061 : IWvithe fc Cenl 120071 ). The observa- 
tion of anomalously large EWs in Lya emitting galaxies at 
high redshift may therefore provide observational evidence 
for such a scenario. In the future, t he He 1640 A may be 
used as a complementary probe (e.g iTumlinson et al.ll200ll . 
12003! ) . The EW of this line is smaller by a factor of ^20 for 
population III j|Schaererl l2003'). However, the He 1640A will 
not be subject to a small transmission of ~ 10 — 30%, mak- 
ing it accessible to the next generation of space telescopes. 
On the other hand, it may also be possible to observe the 
He 1640 A in a composite spectra of z = 6 LBGs. Indeed, 
the He 1640 A line has alread y been observed in t he com- 
pos ite spectrum of z= 3 LBGs l|Shaplev et al.ll2003l ). which 
led [j lmcncz fc Haimanl (|2006l ) to argue for population III 
star formation at redshifts as low as z = 3 — 4. If population 
III star formation was more widespread at higher redshifts, 
as predicted by our model, then the composite spectrum 



their lifetimes. Our model does not specify when this population- 
Ill phase occurs. Hypothetically, the population III phase could 
occur at an arbitrary moment in the galaxies' lifetime when it 
is triggered by a merger of a regular star forming galaxy and a 
dark matter halo containing gas of primordial composition. Note 
however that such a model would probably have difficulties ex- 
plaining the apparent observed correlation between Lya EW and 
the age of a stellar population f § 16.11 1. 
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of LBGs at higher redshifts should exhibit an increasingly 
prominent He 1640 A line. In particular, this line should be 
most prominent in the subset of LBGs that have large EW 
Lya emission lines. 
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